Quantum deflagration and supersonic fronts of tunneling in molecular magnets 
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Theory of magnetic deflagration taking into account dipolar-controlled spin tunneling has been 
applied to the realistic model of molecular magnet Mni2 Ac. At small transverse field, the front 
speed v has tunneling maxima on the bias field B z reflecting those of the molecular spin's relaxation 
rate calculated from the density-matrix equation. At high transverse field, spin tunneling directly 
out of the metastable ground state leads to front speeds that can exceed the speed of sound. Both 
for the weak and strong transverse field, the spatial profile of the deflagration front near tunneling 
resonances shows a front of tunneling that triggers a burning front behind it. 

PACS numbers: 75.50.Xx, 75.45.+j, 76.20.+q 



I. INTRODUCTION 

Burning or deflagration,^] a self-supporting phe- 
nomenon that can exist in the form of propagating fronts, 
is decay of metastable states, controlled by the temper- 
ature increasing as a result of the energy release and 
heat conduction toward the cold region before the front. 
The main ingredient of deflagration is the decay rate 
T of the metastable state that has the Arrhenius form 
T = Tq exp [— U/(kBT)] at low temperatures T <C J7, 
where U is the energy barrier. One could ask if defla- 
gration can exist in magnetic systems, many of which 
are bistable due to a strong uniaxial anisotropy that cre- 
ates an energy barrier between the two energy minima. 
However, the energy release in magnetic systems is much 
weaker than in the case of a regular (chemical) deflagra- 
tion, thus at room temperatures the ensuing tempera- 
ture increase is too small to change the relaxation rate 
and support burning. The situation changes at low tem- 
peratures, however, since temperature generated by de- 
cay of metastable states can exceed the initial temper- 
ature by far and result in a strong increase of T. Re- 
cently magnetic deflagration has been observed in low- 
temperature experiments on the molecular magnet Mni2 
AcP4J This discovery initiated theoretical^ and further 
experimental^ work. Magnetic deflagration has also 
been observed on manganitesP Very fast moving fronts 
of burning in Mni2 Ac initiated by a fast sweep of the 
magnetic field have been observed in Ref. [TUJ This leads 
to the idea of magnetic detona tion d riven by thermal ex- 
pansion creating a shock waveP^HH 

Main exponents of magnetic deflagration, molecular 
magnets, are built of molecules with a large effective 
spin, such as S = 10 in Mni2 and Feg. Their uniaxial 
anisotropy D crea tes the energy barrier DS 2 c± 67K for 
spin rotatiorP^E^ (see Ref. [15] for a review). Molecular 
magnets made quite a big splash by the discovery of res- 
onance spin tunneling 1 -^"^ that occurs when spin energy 
levels on different sides of the barrier match. This is con- 
trolled by the bias created by the longitudinal magnetic 
field. Magnetic molecules in molecular magnets form a 
crystal lattice (body-centered tetragonal for Mni2 Ac). 



As magnetic cores of the molecules are shielded by or- 
ganic ligands, there is no exchange interaction between 
the molecules in the crystal, and the dipole-dipole inter- 
action (DDI) is dominating. Different members of the 
Mni2 family remain in the center of magnetic deflagra- 
tion research because of the elongated shape of the crys- 
tals. To the contrast, Fes crystals have pyramidal shape, 
inappropriate for studying moving fronts. 

The impact of spin tunneling on deflagration in molec- 
ular magnets has been addressed in Refs. 114|5|191 Since 
no transverse magnetic field was applied in experiments 
so far, tunneling via low- lying states was negligibly small. 
Thus quantum effects in deflag ration could only exist due 
to thermally assisted tunneling } 20 * 21 * via the energy levels 
just below the top of the barrier. This effect can be taken 
into account as effective lowering of the barrier U at res- 
onant values of the biasP^ Peaks of the deflagration front 
speed vs longitudinal magnetic field (Fig. 4 of Ref. [4]) 
have been interpreted as spin tunneling. The simplest 
way to explain these peaks was to use the escape rate T 
with the effective barrier U in the standard formula for 



the speed of the deflagration front, Eqs. (10) and (11) 
with v = 1 (dashed line in Fig. 4 of Ref. |4J. For higher 
bias and thinner crystals, observed speed maxima were 
much weaker (Fig. 5 of Ref. 6 and Fig. 3 of Ref. [7j) that 
created a controversy. 

At the same time, there was a quest for an essentially 
quantum mechanism of deflagration in molecular mag- 
nets that does not reduce to mere barrier lowering in 
the thermally activated escape rate. As a further de- 
velopment, fronts of spin tunneling (dubbed "cold defla- 
gration") controlled by the dipolar field at zero temper- 
ature have been proposed! 23 * 24 * This mechanism requires 
a strong transverse magnetic field that creates a suffi- 
ciently large tunnel splitting A between the metastable 
ground state and an excited state on the other side of 
the barrier. The idea is that dipolar field created by the 
sample produces a bias on magnetic molecules (spins) 
that is typically large in comparison to A, thus the dipo- 
lar field can control tunneling. As tunneling of one spin 
changes dipolar fields on other spins, facilitating or pre- 
venting their tunneling, the problem is self-consistent. It 
was shown that there are solutions in which the spatial 
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distribution of magnetization and dipolar field is adjust- 
ing in such a way that there is a moving front of spin 
tunneling with many spins in the front core being on 
resonance that allows them to tunnel efficiently. This 
so-called laminar front has been found for not too large 
values of the external bias. For a larger bias it breaks 
down, resulting in a slow non-laminar front where most 
spins are off-resonance. 24 Fronts of cold deflagration ex- 
ist within the dipolar window of the external bias having 
the width equal to the dipolar field B Z D ^ = 52.6 mT 



duced by a uniformly magnetized molecular magnet ! 25 l 26 l 
In addition to the transverse field, observation of fronts 
of tunneling in pure form requires a good thermal contact 
between the crystal and its environment, so that released 
heat gets conducted away and the temperature remains 
low. 

If the crystal of a molecular magnet is thermally in- 
sulated, spin tunneling in a biased case leads to release 
of Zeeman energy and the temperature increase. In this 
case both spin tunneling and thermal activation can play 
a role, so that deflagration is controlled by two parame- 
ters, dipolar field and the temperature. The combined 
quantum-thermal theory of magnetic deflagration has 
been proposed in Ref. |27l In contrast to the pure cold 
deflagration, where in the case of overdamped tunnel- 
ing it is sufficient to use the Lorentzian form of the tun- 
neling rate near the resonance [Eq. (12) of Ref. [24] . 
here one needs the numerically calculated escape rate 
T(B Z1 T) for both resonant and non-resonant values of 
B z . This escape rate has been calculated from the density 
matrix eq uation ^ based on the universal spin-phonon 
interaction! 29 1 30 1 To the contrast with the pure cold de- 
flagration that leaves some metastable magnetization un- 
burned behind the front, the combined deflagration leads 
to complete burning, as the standard magnetic deflagra- 
tion. This flattens out irregularities of non-laminar fronts 
and makes them move faster, reaching high speeds at the 
right end of the dipolar window (see Fig. 4 of Ref. |27|) . 

Ref. [271 used the generic model of a molecular magnet 
with the anisotropy of the form —DS Z . In this model, 
tunneling resonances of all levels take place at the same 
value of B z 

B z = B k = kD/( gl jLB), k = 0, ±1, ±2, . . . (1) 

and, non-trivially, the resonances remain unchanged if 
transverse magnetic field is applied. In the real Mni2 Ac 
there is an additional term —AS^ that makes resonances 
of different levels be achieved at different values of B z . 
The latter was used to experimentally monitor the tran- 
sition between th erma lly assisted and ground-state tun- 
neling in Mni2 Ac ! 31 * 32 * Splitting of tunneling resonances 
should manifest itself in experiments on magnetic defla- 
gration, and studying related phenomena is one of the 
aims of this work. 

Another aim of this work is to explore the high-speed 
regime of magnetic burning near the ground-state reso- 
nance at high transverse fields. As the speed of fronts of 
tunneling should be much higher than that of the stan- 



dard burning fronts driven by heat conduction, burning 
in these fronts should be independent of the thermal dif- 
fusivity, that resembles detonation. To study this regime, 
more accurate numerical calculations on longer crystals 
have to be performed. 

The rest of the paper is organized as follows. In Sec. [TT] 
equations describing deflagration with dipolar-controlled 
spin tunneling are set up and the method of their solution 
is outlined. Sec. [Till introduces the relaxation rate of 
magnetic molecules that is calculated with the help of the 
density-matrix formalism and contains the effects of both 



thermal activation and spin tunneling. Sec. IV presents 
numerical results for the front speed in weak transverse 
fields. Sec. |V|is devoted to the case of a strong transverse 
field, where ground-state tunneling leads to supersonic 
front speeds. Concluding section summarizes the results 
obtained and outlines unsolved problems. 



II. EQUATIONS OF DEFLAGRATION WITH 
SPIN TUNNELING AND DIPOLAR FIELD 



The system of equation describing deflagration with 
quantum effects in molecular magnets 27 consists of the 
rate equation for the metastable population n 



dn(t, z) 
dt 



T(B tot , z (z),T(z)) n(t,z)-n^(T) (2) 



and the heat conduction equation that can be conve- 
niently written for the thermal energy £ per magnetic 
molecule 



d£(t,z) 
dt 



d_ d£(t,z) 
dz dz 



A£^M. (3) 



It is assumed that the crystal has an elongated shape 
and everything depends only on the coordinate z along 
the geometrical axis of the crystal. The easy axes of mag- 
netic molecules are also directed along this axis, that was 
the case for all experimentally studied crystals E£I4I£M|26| 
In Eq. Q Y{B Z ,T) is the numerically computed relax- 
ation (escape) rate of magnetic molecules' spins out of 
the metastable state with the spin pointed to the left, 
when a longitudinal external field is applied in the di- 
rection to the right. n^ eq \T) is the thermal-equilibrium 
population of the metastable state that is small in the 
case of a large bias and will be discarded. In Eq. ([3| k is 
thermal diffusivity that proves to be difficult to measure. 
Estimations 3 yield k ~ 10 _5 m 2 /s (comparable to that of 
metals) that will be adopted here. The second term in 
this equation is the source term, in which AE is the en- 
ergy released by transition of one molecular spin from the 
metastable state to the ground state |— S) —> \S), that 
is, AE — 2SgjiBB z . The relation between the energy £ 
and temperature is given by 



£{T) 



f 

Jo 



C(T')dT', 



(4) 
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where C(T) is the experimentally measured heat capacity 
of Mni2Ac per magnetic moleculeP^ 

Since the relaxation rate T(B totjZ ,T) has very sharp 
maxima at the resonance values of the total longitudinal 
field B totiZ , it is important to include the dipolar field 
created by the crystal, 



D 



B tot)Z (z)=B z + B { z D \z). 



(5) 



Although the dipolar field B Z D ^ is much weaker than the 
external field B z (and thus can be dropped in AE), it 
is much greater than the width of tunneling peaks in 
r(B to t,ziT), so that it can control tunneling. It is con- 
venient to represent B Z D ^ in the form 



Sg^i 



D z 



(6) 



where D zz is the dimensionless dipolar field, = 
is the unit-cell volume, a and c are lattice spacings. 
For Mni2 Ac one has Sg/iiB/vo = 5.0 mT. For crys- 
tals of cylindrical shape with radius R and length L one 
obtain^ 



D zz (z 



r±j c 

) = / dz' 



2iruR 2 a z (z / ) 



(z' - z) z + R 2 



3/2 



k D a z (z), (7) 



where v is the number of molecules per unit cell, v = 2 
for Mni2 Ac, o z = 1 — 2n is polarization of pseudospins 
representing spins of magnetic molecules (a z = ±1 in the 
ground and metastable states, respectively) and 
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Sttu/3 - = 4tt^ - > 0, 



(8) 



Here the barred quantities correspond to the reduced 
dipolar field inside a uniformly magnetized sphere and 
a long cylinder, and D zz = D zz a z for a z = const. For 
Mn i2 Ac calculations yielcP L>i? h) = 2. 155, D z 7 l) = 
10.53 (in real units B Z D) = 52.6 mT™), and thus in 
the local term of Eq. ^ one has k& = 14.6. One 
can check that Eq. ^ yields the correct result for 
the field inside a long uniformly magnetized cylinder. 
At the ends of a cylinder dipolar field has the form 

D zz = (j} z s z h ^ — 27r^/3^ a z that for Mni2 Ac becomes 

D zz = — 2. 03(72 . Dipolar field opposite to the spin 
orientation is the reason for the instability of the uni- 
formly magnetized state of Mni2 Ac that leads to domain 
formation. 34 For other shapes such as elongated rectan- 
gular, one obtains qualitatively similar expressions.^ 

It has to be stressed that the results above represent 
the dipolar field exactly at the magnetic molecules in the 
lattice and they depend on the lattice structure. Using 
the spatially averaged field following from macroscopic 
magnetostatics would be a mistake. Indeed, the magne- 
tostatic field inside a long uniformly magnetized cylinder 
is B Z D ^ = 47rM, where the magnetization is given by 



M = vSqiib/vq. The dipolar field in Mni2 Ac is es- 
sentially smaller, B Z D) = (D z c z yl) /2)M = 5.26M. The 
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Figure 1: Reduced dipolar field in a deflagration front in 
the slow-burning limit, created by the magnetization profile 
a z (z) = — tanh [(z — zo)/ld]. Solid line: Eq. Q; Points: Di- 
rect summation of dipolar fields over the Mni2 Ac lattice. 



difference between the two is due to the local term with 
k D in Eq. Q. 

A striking feature of Eq. ([7]) is that the integral and lo- 
cal terms have different signs. The integral term changes 
at the scale of R while the local term can change faster, 
that creates a non-monotonic dependence of D zz (z). In 
the case of a regular magnetic deflagration, the spa- 
tial magnetization profile in the slow-burning limit is 
&z( z ) = — tanh \{z — zo)/ld\, where Id is the width of the 
deflagration front that satisfies Id <C R (see below). The 
resulting dipolar field is shown in Fig. [T] where the line 
is the result of Eq. ^ and points represent the dipolar 
field along the symmetry axis of a long cylindrical crys- 
tal calculated by direct summation of microscopic dipolar 
fields over the Mni 2 Ac lattice. One can see that Eq. ^ 
is pretty accurate, small discrepancies resulting from Id 
being not large enough in comparizon to the lattice spac- 
ing a. The central region with the large positive slope 
is dominated by the local term of Eq. ^ that changes 
in the direction opposite to that of the magnetization. 
For R ^> Id, D zz reaches the values ±14.6 due to the 
local term before it begins to slowly change in the oppo- 
site direction. In real units the dipolar field at the local 
maximum and minimum is ±B z kD \ where 

£f D) =72.9mT (9) 

exceeding the dipolar field of the uniformly magnetized 
long cylinder B Z D ^ = 52.6 mT. Also one can see from 
Fig. [7] that the dipolar field becomes opposite to the 
magnetization at the ends of the cylinder, as mentioned 
above. 

Equations (|2|-([7|) form a system of integro-diflerential 
equations describing deflagration with spin tunneling in 
molecular magnets taking into account the dipole-dipole 
interaction. Before discussing numerical solution of these 
equations, it is worth recuperating the results of the stan- 



dard ("hot") deflagration and of the cold deflagration. If 
the whole released energy remains in the body and the 
initial temperature is very low, the thermal energy per 
spin behind the front is AE. The corresponding temper- 
ature defined by the inversion of Eq. Q is the so-called 
flame temperature Tf = T(AE) that is in the range 10- 
15 K in deflagration experiments. Theory of deflagration 
yields the expressions for the speed of the front v and 
front width Id 



v = vr f l d = vy/KfTf, ld = yj «//r/, (10) 

where and 

T f = T exp (-W f ) , W f = U/(k B T f ) (11) 

are thermal diffusivity and relaxation rate at the flame 
temperature, while v is a dimensionless coeficient. It was 
showrP that in the slow-burning limit Wf ^> 1 one has 
v = 2/ \/Wf. On the other hand, the speed of the laminar 
front of tunneling at zero temperature is given b y 23 ! 24 ! 

v = v*T res R, (12) 

where r res = A 2 /{h 2 T m >) is the relaxation rate at over- 
damped tunneling resonance, A/h< T m /, T m > being the 
decay rate of the matching level m! at the other side 
of the barrier, R is the width of the crystal (radius of 
the cylinder in our model), and v* is a dimensionless 
coefficient. With a sufficiently strong transverse field ap- 
plied, one can have A/h ~ T m t at the applicability limit 
of the overdamped approximation, and then r res ^> Tf 
because thermal activation goes over high levels of the 
magnetic molecule where the distances between the lev- 
els and thus the energies of phonons involved are much 
smaller than for the low-lying levels, and also because 
is exponentially small since Tf < U. Additionally, 
estimation of Id with Kf = 10 _5 m 2 /s and the experi- 
mental value To = 10 T s _1 yield Id ~ 3 x 10 -4 mm for 
B z near the first tunneling resonance and even smaller 
for larger bias. As in the experiment the width of the 
crystal was much larger than Id (0.3 mm in Ref. El 0.2 
mm in Ref. [6] and 1 mm in Ref. 4), one can see that 
T res R ^> T fid is quite possible in a strong transverse 
field, and then the front of spin tunneling is much faster 
than the front of spin burning. A very conservative es- 
timation with r res To = 10 7 s _1 and v* 1 for the 
crystal 0.2 mm thick yields v ~ 1000 m/s. As said above, 
in a strong transverse field one can have r res ^> To, so 
that the speed of a spin-tunneling front can easily surpass 
the speed of sound that is about 2000 m/s in molecular 
magnets (see analysis in Ref. [35j). The results of our 
calculations confirm this. 

Discretization of the variable z reduces Eqs. [2j [3j and 
[7] to a system of ordinary differential equations that can 
be solved numerically. Very narrow tunneling peaks in 
T(B Z1 T) make it necessary to carefully control the step 
in the numerical integration. Wolfram Mathematica's 
NDSolve proves to be an efficient tool for this problem. 
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Figure 2: Relaxation rate of Mn^Ac vs temperature and lon- 
gitudinal magnetic field in the transverse field B± = 0.04 T. 
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Figure 3: Relaxation rate of Mn^Ac vs temperature and lon- 
gitudinal magnetic field in the transverse field B± = 0.04 T 
for a stronger bias. 



To ignite a deflagration front, the temperature at the 
left end of the crystal had been increased during a short 
time. Then the equations were solved and, to find the 
front speed, the time of arrival of the front at the right 
end of the crystal was measured. 



III. THE RELAXATION RATE 

It is crucial to calculate and tabulate the relaxation 
rate T(B Z1 T) before solving the deflagration problem be- 
cause a runtime calculation of T(B Z , T) is practically im- 
possible. We use the effective-spin model with the Hamil- 
tonian containing the uniaxial anisotropy —B)S 2 Z — AS Z 
and other anisotropy terms, according to Ref. |36[ Spin- 
phonon interaction is taken into account within the uni- 
versal model of pure rotations of the crystal field by 
transverse phonons described in Refs. 29 30 37. Since 




Figure 4: Relaxation rate of Mn^Ac vs temperature and Ion- Figure 5: Magnification of the ground-state tunneling peak of 
gitudinal magnetic field in the transverse field B± = 3.5 T. F(B Z ,T) (multiplied by 100) at B± = 3.5 T. 



in this model the crystal field is not distorted, spin- 
phonon coupling coefficients can be expressed through 
the measurable crystal- field parameters. The density- 
matrix equation has been solved within the semi-secular 
approximation that is valid everywhere, including tun- 
neling resonances! 2 ^ 

In the generic model of a molecular magnet with the 
anisotropy —DS^ the fields corresponding to tunneling 
resonances are given by Eq. for all level pairs. The 
resulting T(B Z ,T) in a strong transverse field is shown 
in Fig. 2 of Ref. |27l Tabulation of such a function re- 
quires a lot of points along the B z axis in the vicinity of 
tunneling maxima. The realistic model with the uniaxial 
anisotropy —DS 2 Z —AS\ is more complicated because tun- 
neling resonances for different level pairs are achieved for 
different B z that depend on the transverse field. Thus the 
first step is to find tunneling peaks numerically for a given 
transverse field, then to build a non-equidistant grid with 
a small step near the peaks, then calculate T(B Z ,T) and, 
finally, make the interpolation. These tasks have been 
fulfilled with the help of Wolfram Mathematica using a 
high custom precision and parallelization. 

For a weak transverse field (set to B± = 0.04 — 0.05 
T that may result from a 1° misalignment between the 
crystal axis and the longitudinal field) T(B Z ,T) contains 
a zoo of tunneling peaks shown in Fig. [2] The range 
of B z here corresponds to that in Ref. [4 and contains 
groups of resonances with k = 2, 3 and partially 4. One 
can see that ground-state resonances, that are the only 
survivors at T = 0, are achieved at higher fields than 
resonances of excited states. At temperatures as high 
as flame temperature, low-lying tunneling resonances are 
drowned in the non-resonant background. There is also 
a much weaker non-resonant tunneling at T = 0. Relax- 
ation rate at a stronger bias, also in a small transverse 
field, corresponding to that in Refs. I3)6f f8| is shown in 
Fig. [3] At such bias, the effect of ground-state tunneling 
begins to appear at high temperatures. 



In a strong transverse field such as B± = 3.5 T in Fig. 
[4] the barrier is strongly lowered and most of tunneling 
resonances are broadened away. Here one can see the 
ground-state resonance (B z = 0.522), the first-excited- 
state resonance (B z = 0.490), and with an effort a very 
broad second-excited-state resonance further to the left. 
Note the much higher tunneling rate at T = 0, in com- 
parison with the previous figure. The range of B z in Fig. 
[4] corresponds to the tunneling resonance with k = 1. 
Fig. [5] shows the details of the ground-state peak in Fig. 
[4] The height and width of this peak increase with tem- 
perature. This increase is moderate, however, in com- 
parison to the exponential increase of the non-resonant 
relaxation rate. The first-excited-state peak in Fig. [4] 
is higher than the ground-state peak at the flame tem- 
perature but is plays a much smaller role in the front 
propagation, as we will see below. 

A long-standing problem in the theory of relaxation 
of molecular magnets is the prefactor To in the Arrhe- 
nius relaxation rate being by two orders of magnitude 
too small. This was already recognized in the early Ref. 
21, Without introducing artificially strong spin-phonon 
interactions it is imp ossible to arrive at r ~ 10 7 s _1 
observed in experiment d 22 * 33 * using the standard spin- 
lattice relaxation model considering one spin in an in- 
finite elastic matrix. This model could be justified for 
a strongly diluted molecular magnet but in the normal 
case it cannot. High density of magnetic molecules should 
lead to such collective effects as superradiance^EU an d 
phonon bottleneck! 42 " 44 ! As it would be difficult to deal 
with these complicated issues while addressing the quan- 
tum deflagration problem, the calculated relaxation rate 
was simply multiplied by 100 to approximately match 
the experiment. It is instructive to plot the theoreti- 
cal deflagration speed given by Eq. ( [To| ) (with v = 1) 
at small transverse field as function of B z using the cor- 
rected values of T(B z ,Tf(B z )). Fig. [6] shows a good over- 
all agreement, except for tunneling maxima in the micro- 




Figure 6: Speed of deflagration front estimated from 
Eq. (10) for the microscopically calculated relaxation rate 
T(B Z , Tf (B z )) with the correction factor 100, together with 
Eq. (flop using the Arrhenius formula for T(B z ,T f (B z )). 



scopically calculated result. As tunneling resonances are 
broadened by ligand disorder, dipolar field, and nuclear 
spins, very narrow peaks due to tunneling resonances of 
lower levels here will be washed out in the experimentally 
measured front speed. In fact, a similar interpretation of 
experimental results have been done in Ref. 3J where 
the dashed line in Fig. 4 is ^KfY f with Tf taken from 
relaxation experiments on the same crystal. 

An alternative explanation of much higher relaxation 
rates observed in the experiment is based on deviations 
from the strong-exchange model that lead to mixing of 
the states with different total spin S. In Ref. [45] it was 
shown that this small mixing taken into account pertur- 
batively leads again to the giant-spin model with S = 10, 
however, with additional higher-order crystal- field terms 
that would normally be absent for <i-electrons. These 
additional terms can explain the observed ground-state 
tunnel splitting A in Feg that is three orders of mag- 
nitude larger than the theoretical result using the stan- 
dard spin Hamiltonian. Similar mechanism could work 
for Mni2 and lead to the increase of the spin-lattice re- 
laxation rate as well. However, the importance of this 
mechanism is limited to small transverse fields. Most 
interesting results below for supersonic fronts of tunnel- 
ing directly out of the metastable ground state without 
thermal activation require a strong transverse field that 
produces a large tunnel splitting. In this limit, the latter 
becomes insensitive to crystal-field terms responsible for 
tunneling in zero or small transverse fields. 



IV. FRONT SPEED AT WEAK TRANSVERSE 
FIELD 

The procedure of numerical solution of the quantum 
deflagration equations is discussed at the end of Sec. [I[ 



Figure 7: Numerically calculated speed of the deflagration 
front in a Mni2 Ac crystal in small transverse field. 
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Figure 8: Numerically calculated speed of the deflagration 
front for small transverse field for a Mni2 Ac crystal with a 
smaller transverse size. 



The result for the front speed at small transverse fields 
in the range B z =0.7-1.7 T is shown in Fig. [F| Here 
the cylinder radius R in our model has been chosen so 
that it yields the same cross-section as the crystal of 
transverse sizes L a = = 1 mm in Ref. [4j that is, 
R = ^LaLb/it = 0.564 mm. One can see that, in com- 
parison to Fig. |6j narrow tunneling peaks are washed 
out and only broad peaks remain. The reason is that 
the total magnetic field in the crystal is not constant and 
changes in the front as shown in Fig. [T] so that tunnel- 
ing resonances in v are spread. Overall there is a good 
agreement between our Fig. [2] and Fig. 4 of Ref. |4j For 
a comparizon, the calculated front speed for a crystal of 
smaller transverse dimensions, L a = = 0.2 mm, such 
as in Refs. [BHHlis shown in Fig. [8] In this case tunneling 
peaks are not washed out, although they are much wider 
and lower than those in Fig. [6] Some of these peaks are 
asymmetric, similarly to the single large peak in Fig. 4 
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Figure 9: Numerically calculated speed of the deflagration 
front in a thinner Mni2 Ac crystal for small transverse field 
and a larger bias. 



of Ref. [27| The reason for this asymmetry will be dis- 
cussed below. Then, Fig. [9] shows the calculated front 
speed for the bias and crystal size corresponding to the 
experiments in Refs. [6HSJ Here tunneling peaks are quite 
pronounced, at variance with the above experiments that 
show very small peaks. Just above 3 T and just below 
3.5 T there are regions where the speed is too high to 
be measured in this calculation, an effect of ground- state 
tunneling. 

Spatial profiles of the magnetization, energy, and the 
total bias field in the deflagration front give an idea of 



the role played by spin tunneling. Fig. [10] shows these 
profiles at B z — 1.5 T that is far from resonances. In this 
case there is a pure slow burning with the magnetization 
and energy profiles of a tanh shape. 5 Dipolar field shown 
in the lower panel plays no role in the process. 

Fig. 11 shows the spatial profiles at the asymmetric 
peak of v at B z = 2.852 T in Fig. (9) Here the front 
speed is high because of tunneling at the face of the front 
where in the lower panel the total bias field is flat at 
the level of the tunneling resonance at B z ^ tot = 2.889 T. 
Magnetization distribution adjusts so that the dipolar 
field ensures resonance for a sizable group of spins that 
tunnel. Tunneling of these spins results in energy release, 
the temperature and relaxation rate increase, and tunnel- 
ing gives way to burning in the central and rear areas of 
the front. 

Formation of the asymmetric maxima of the front 
speed can be explained as follows. When B z increases, 
the peak of B z ^ tot that arizes due to the local dipolar field 
reaches the resonant value. Here the strong increase of 
v(B z ) begins. The maximum of B Zjtot sticks to the reso- 
nance value and becomes flat with progressively increas- 
ing width. Greater width of the resonance region results 
in a stronger tunneling and higher front speed. With fur- 
ther increase of B z , the right edge of the tunneling region 
moves too far away from the front core into the region 



where the temperature is too low. As the tunneling res- 
onance in question is thermally assisted, it disappears at 
low temperatures, thus the flat region of B z ^ tot cannot 
spread too far to the right. As a result, the flat config- 
uration of B z ^ tot becomes unstable and suddenly B z ^ tot 
changes to the regular shape of Fig. [T] that crosses the 
resonance twice in the face part of the front. At the right 
crossing the temperature is too low and tunneling does 
not occur, whereas at the left crossing burning already is 
going on and tunneling cannot add much. There can be 
the third resonance crossing further to the left but it does 
not play a role because everything has already burned. 
It should be noted that multiple resonance crossings do 
not occur in the laminar regime of the pure quantum case 
(cold deflagration), see Fig. 2 of Ref. [24j 

If the transverse size of the crystal is large, R ^> 
the slope of B Zjtot to the right of the maximum in Fig. W\ 
is small. In this case increasing B z leads to a very quick 
displacement of the right border of the tunneling region 
to the right where tunneling cannot take place, as ex- 
plained above. Thus tunneling peaks of v(B z ) should be 
very narrow for such crystals. This explains why tunnel- 
ing peaks are quite pronounced in Fig. [8] but very small 
in Fig. [7| 

Tunneling peaks of v(B z ) corresponding to broad res- 
onances of highly excited states are almost symmetric, 
such as the high peak at B z = 2.644 T in Fig. [9] In 
this case peaks are formed when the maximum of B Zjtot 
crosses the resonance. In these cases progressive flat- 
tening of B Zytot does not occur because here tunneling 
requires high temperatures and the right border of the 
tunneling region cannot move to the cold region to the 
right. 

Fig. [T2| shows an off-resonance front again, B z = 3 T, 
but it is not a slow burning front anymore because the 
bias is high. In this region the analytical theory of Ref. 
5 does not work that can be seen on the magnetization 
and energy profiles that differ from the tanh shape. 



V. FRONT SPEED AT STRONG TRANSVERSE 
FIELD 

As one can see from Fig. [4j at strong transverse fields 
the structure of the relaxation rate T(B Z ,T) simplifies 
because tunneling resonances of the most excited states 
broaden away. At B± = 3.5 T one can see only two tun- 
neling peaks, and the ground- state tunneling peak is not 
drowned by the thermal-activation processes up to the 
highest temperatures. This means that in a bias win- 
dow around this peak, the barrier is cut completely. The 
latter changes the dynamics of the system, drastically 
increasing the role of tunneling in the front propagation. 
Since tunneling out of the metastable ground state does 
not require an elevated temperature, the right border of 
the tunneling region before the main part of the front 
can shift unlimitedly to the right without causing the in- 
stability that kills tunneling, described in the preceding 
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Figure 10: Spatial profiles of the deflagration front in a small 
transverse field, B± — 0.05 T. The bias B z — 1.5T is far from 
resonances, thus the front is that of a pure slow burning. 



Figure 12: Spatial profiles of the deflagration front in a small 
transverse field, B± — 0.05 T. Faster burning with no tunnel- 
ing. 
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Figure 11: Spatial profiles of the deflagration front in a small 
transverse field, B± = 0.05 T at the peak of the front speed 
at B z = 2.852 T. There is a resonance spin tunneling at the 
face of the front and burning in its central and rear parts. 



Figure 13: Quantum self-ignition before a slow-burning front, 
leading to a fast combined tunneling-burning front near 
ground-state resonance at high transverse fields. 



section. Thus the width of the tunneling region can reach 
the values of order i? p 3 l 24 l that leads to front speeds much 
greater than the speed of a regular magnetic deflagration 
[see comments after Eq. (12)]. 

Numerical results at high transverse fields show that 
shortly after ignition by raising the temperature at the 
left end of the crystal, a regular slow-burning front 
can transform into a fast combined tunneling-burning 
front by quantum self-ignition before the slow-burning 
front, if the crystal is near ground- state tunneling reso- 
nance. Fig. [13] shows this phenomenon at B± = 3.5 T 
and B z = 0.47 T, where the ground-state resonance is 
achieved at B z ^ tot = 0.522 T. One can see that at short 
times, tT f ~ 10, there is a slow front with a steep profile 
but before the front, where B Zjtot crosses the resonance 
value, spins begin to tunnel. This quantum self-ignition 



leads to flattening of the B Z}tot curve and formation of an- 
other, fast moving front, with tunneling followed by burn- 
ing. The spatial profile of B Zytot at different times that 
shows that self-ignition before the slow-burning front is 
caused by spin tunneling, is shown in Fig. [14] 

The front speed v in the vicinity of a biased ground- 
state resonance in a strong transverse field can achieve 
supersonic values, as can be seen in Fig. 
in accord with the comments below Eq. 



15 



(12). 



This is 
To the 



contrast, the small peak on the left side of Fig. 15 is due 
to the first-excited-state tunneling resonance at Bk = 
0.490, see Fig. [4] Its position is given by B z — B^ — 
g( k n) _ o 417 x that is close to the position in the figure. 
In Sec. IV it was explaind that front of tunneling via 
excited levels cannot shift much ahead of the burning 
zone because it is too cold before the front. This limits 
the speed of such fronts and explains why the speed of the 
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Figure 14: Profiles of the total longitudinal magnetic field at 
B± — 3.5 T and B z — 0.47 T at different moments of time. 
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Figure 15: Front speed for a strong transverse field (B± =3.5 
T) in the vicinity of the ground-state tunneling resonance at 
0.522 T. There is a strong increase of the front speed within 
the dipolar window of 125.5 mT around the resonance. The 
small peak on the left (inset) is due to the first-excited-state 
tunneling resonance. 



first-excited state tunneling front is much smaller than 
that of the ground-state tunneling front, in spite of the 
relaxation rate at the former being higher. 

Returning to the ground-state tunneling front, it 
should be stressed that no metastable population is left 
behind the front (see Fig. [l3] ), although there is unburned 
metastable population behind pure non-thermal fronts of 
tunneling. 23 24 Here, the metastable population is burn- 
ing just behind the front of tunneling as the result of 
the temperature increase. It should be stressed that heat 
conduction cannot support burning fronts moving faster 
than the speed of sound and it becomes non-operative 
in this case. In this respect, the situation is remind- 
ing of detonation that has been suggested for molecular 
magnets in Refs. lll|12l in the case of a strong bias and 
thus high energy release. In detonation, thermal expan- 



sion resulting from burning sends a shock wave into the 
cold region before the front where, as a consequence, the 
temperature rises as a result of compression, initiating 
burning. As the mechanism of detonation is based on 
elasticity, the speed of a detonation front is comparable 
to the speed of sound. Fronts of tunneling are not based 
on elasticity and their speed can be much higher. How- 
ever, shock waves must accompany tunneling fronts and 
modify their properties in some way. Experimentally, fast 
deflagration or detonation fronts in Mni2 Ac have been 
observed in Ref. [10] but they were caused by a very fast 
sweep of B z , so that there is a question to which extent 
the process was self-propelled. 
One can see in Fig. 



15 that the speed of the front 



is asymmetric and grows toward the right end of the 
tunneling window, showing divergence or nearly diver- 
gence of the front speed. In the case of cold deflagra- 
tion (assuming the unbroken laminar regime everywhere) 
v diverges at the right border of the dipolar window, 
B z — Bk + B Z D ^ , where Bk is the resonance field of Eq. 
Jli. It is given 



v ~ T re .R- 



B z — Bh 



Bh 



for Bk < B z < Bk + B Z D \ whereas above Bk + B 
it abruptly drops to a zero. The reason for this is that 
above Bk + B Z D \ the total field well before the front 
is above the resonance, so that resonance crossing can- 
not occur. To the contrast, just below Bk + B Z D ^ the 
field well before the front is a little bit below the res- 
onance and increases closer to the front. In this case, 
there is a wide region where the system is close to the 
resonance, and the front speed becomes very high. Thus 
as B z crosses the value Bk + B Z D ^ from below, the front 
speed drops abruptly. Similar behavior can be seen in 



B z 



(13) 



AD) 



Fig. 15 The ground-state resonance is at 0.522 T, and 



adding the dipolar field B 



0.0526 T one obtains 



0.573 T, as in the figure. However, here v drops to the 
speed of the regular magnetic deflagration, as also in Fig. 
4 of Ref. [27j Another difference is that in the case of cold 
deflagration^ tunneling begins at B z = Bk (left border 
of the dipolar window) whereas in our case it begins when 
the local maximum of B Zjtot first touches the resonance 



(see, e.g., the lower panel of Fig. 11). Since for the crys- 
tals studied here the front width is much smaller than 
the transverse size of the crystal, the dipolar field at the 
maximum is close to B z k °^ given by Eq. Thus the 

left border of the dipolar window is at B z = Bk — B z Ud \ 
in Fig. [I5lat B z = 0.45 T. The total width of the dipolar 
window of the ground-state tunneling resonance in Mni2 
Ac is 



AB^ D) = B^ D) + Bf d) = 125.5 mT 



(14) 



that is much greater than dipolar windows of excited- 
state tunneling resonances, see, e.g., Fig. [9] 
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In the case of cold deflagration, there is an unburned 
metastable population in the final state behind the front, 
Eq. (41) of Ref. [22] that can be rewritten as 



n f 



B z - B k 



B 



(15) 



(n = 1 before the front). One can see that the change of 
n across the front An = 1 — rif goes to zero at the right 

border of the dipolar window Bk + B Z D ^ . This reconciles 
the situation with the general requirement that the rate 
of change of the magnetization of the crystal M, limited 
by the tunneling parameter A, remains finite. Indeed, 



M oc (1 - n f )v T res R 



B z - B k 



B 



(16) 



reaches only a finite value M oc T res R at the right border 
of the dipolar window before it drops to zero. In the 
present case of tunneling followed by complete burning, 
M is not limited by A and can achieve very high values 
at the right border of the dipolar window. 



VI. DISCUSSION 

Numerical calculations for deflagration fronts with 
dipolar-controlled spin tunneling for the realistic model of 
Mni2 Ac performed in this work have shown many quan- 
tum peaks in the dependence of the front speed v on the 
external magnetic field B Z1 if a zero or small transverse 
field is applied. The multitude of peaks results from the 
splitting of the tunneling resonance by the —AS Z term 
in the crystal field of the magnetic molecule, and peaks 
in v(B z ) reflect those in the relaxation rate T{B Z ,T) of 
the metastable states of Mni2 Ac molecules. The peaks 
of v(B z ) are more pronounced for crystals of a smaller 
transverse size. 

Whereas the results of the calculations for thicker crys- 
tals in the range of smaller bias are in a qualitative accord 
with the experiments of Ref. 4, the results for thinner 
crystals and stronger bias show much stronger tunneling 
peaks in v(B z ) than it was observed in Refs. 6-8. One 
can try to explain the lack of peaks in the experiment 
by the spread of tunneling resonances as the r esult of 
ligand disorder that is pretty strong in Mni2 Ac ! 46 * 47 * It 



has been shown that static disorder that is weaker than 
the dipole-dipole interaction does not destroy fronts of 
tunneling since the magnetization distribution can ad- 
just so that many spins in the front core are still on res- 
onance and can tunnel. 23 However, static disorder that 
is stronger than the DDI cannot be accomodated by the 
latter and should result in spread and suppression of tun- 
neling maxima in v(B z ). The best way to deal with this 
problem is to make experiments on the members of Mni2 
family that do not have ligand disorder. 

Possibly there is a more fundamental reason for the 
near absence of tunneling peaks in the experiments of 
Refs. 6-8 [ The prefactor To in the theoretical relaxation 
rate being by a factor 10 2 smaller than the measured pref- 
actor suggests collective relaxation processes such as su- 
perradiance and phonon bottleneck that can be expected 
in a dense magnetic system such as molecular magnet. 
Collective boosting of relaxation processes should not af- 
fect tunneling, however. Thus the non-resonant back- 
ground in T(B Z1 T) can move up by a factor 100, partially 
drowning tunneling peaks. On the other hand, in our cal- 
culation we have boosted the whole function T{B Z ,T), 
including tunneling peaks. 

Whereas it is impossible to develop a collective the- 
ory of relaxation in molecular magnets in this paper, one 
can modify our density-matrix calculation by multiplying 
spin-phonon coupling amplitudes by 10. This should re- 
sult in the increase of the non-resonant part of T(B Z ,T) 
by 100, while one could expect tunneling peaks to be 
much less changed. Such a calculation has been per- 
formed but its results do not show a strong suppression 
of tunneling peaks in v(B z ). The likely reason for this 
is that at small transverse field tunneling peaks are due 
to thermally assisted tunneling that also gets boosted by 
an artificial increase of spin-phonon interactions. 

Calculations in the case of a strong transverse field, 
making tunneling directly out of the metastable ground 
state operative, show an increase of the front speed 
within the tunneling window around the tunneling reso- 
nance up to supersonic values. It would be highly inter- 
esting to perform experiments on deflagration fronts in 
this region. 

This work has been supported by the NSF under Grant 
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